Localization in an external electric field 
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The impact of an electric field on the electron localization problem is studied within the frame- 
work of a field-theoretic formulation. The investigation shows that the impact of the electric field 
on the localization corrections is governed by the interplay between two time scales, one set by the 
electric field, and the other by the phase relaxation rate. At very low temperatures the scaling 
of the conductivity is governed by the electric field. In this regime the conductivity depends log- 
arithmically on the field, and an arbitrarily small electric field delocalizes the electron states. At 
higher temperatures the behavior of the conductivity is governed by the temperature scaling. In 
this regime the field has no impact on the observable leading localization corrections. 
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I. INTRODUCTION 

The physics of weakly disordered electron systems has 
been the subject of considerable theoretical and experi- 
mental interest over the past years. According to semi- 
classical physics, the electrons in such systems move bal- 
listically in between occasional scattering events due to 
impurities or phonons, which results in diffusive dynam- 
ics. The semiclassical picture has long been known to 
be incorrect in one-dimension, where arbitrarily weak 
disorder leads to all electronic states being localizedi 
A completely new understanding of transport in two- 
dimensions was obtained in the late 1970s by the scaling 
theory of Abrahams et ali^ These authors showed that 
quantum interferences lead to the localization of all elec- 
tronic states in two-dimensions as well, no matter how 
weak the disorder. Accordingly, all two-dimensional sys- 
tems are, strictly speaking, insulators.'^ 

The same conclusion was reached by a field-theoretical 
approach that was pioneered by Wegner.^ Using the 
replica trick to deal with the quenched disorder, he de- 
rived a generalized nonlinear cr-model whose coupling 
constant is the electrical conductivity. Wegner's orig- 
inal theory used bosonic fields, but later applications 
of his method derived the same model starting from a 
fermionic formulation.^ A supersymmetric formulation, 
which avoids the replica trick, was also given,-i 

A third approach to the localization problem were 
mode- mode coupling theories. While the original for- 
mulation by Gotze* missed the interference effects that 
lead to localization in two-dimensions, a self-consistent 
diagrammatic approach by VoUhardt and Wolfle,^ as well 
as an improvement on Gotze's original method, yielded 
results in agreement with both the scaling theory and the 
field-theoretic approaches. They also agree well with nu- 
merical simulations and with experiments. 

All of the above approaches focus on the electronic 
diffusion coefficient, and make a connection to the elec- 



trical conductivity by means of an Einstein relation. By 
contrast, experimentally the conductivity is measured by 
applying an electric field. This raises the question how 
an externally applied electric field affects the localization 
phenomena. 

In the literature, there is no clear answer to this ques- 
tion. Exact results are available for one-dimensional sys- 
tems onlyjiiiiSiii These are characterized by a critical 
electric field Fc such that for field strengths F < Fc the 
states remain (power-law) localized, while for F > Fc the 
states are extended. 

Unfortunately, these exact calculations for one- 
dimensional systems do not reveal how the electric field 
affects the quantum interference effects. The impact of 
an electric field on the localization in two-dimensional 
systems is therefore much less clear. It has been argued 
that a weak homogeneous electric field has no influence 
at all on the Iocalizationpi^*iSii& Other authors have con- 
cluded that an arbitrary small electric field already leads 
to delocalizationiiiiS In between these two extremes lies 
a theory that predicts a strong modification of the weak- 
localization corrections by an electric fieldfiSiSli and the- 
ories in which delocalization only occurs if the electric 
field exceeds a critical value iSLSSiS^. 

Experimentally, the situation is not clear either. While 
early experiments^iiSi seemed to find modifications of the 
weak-localization corrections by a field, no impact of the 
electric field on the localization corrections was observed 
in Ref. Similarly, Ref. concluded that an electric 
field has no impact on the localization corrections. On 
the other hand, the scaling argument of Ref. Ill, which 
implies that an arbitrary small electric field leads to de- 
localization, is often used in the interpretation of exper- 
iments, as discussed in Ref. l28l 

In the present paper we revisit this problem by means 
of field-theoretic techniques. We conclude that delo- 
calization does indeed occur for arbitrarily small fields. 
However, this effect is not observable at realistic tempera- 
tures, which explains some of the apparent contradictions 
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between theory and experiment. We also find that the 
leading contribution to the electric field scaling, which 
was assumed in Ref. to be the dominant efi^ect, has a 
zero prefactor, and the actual effect of the electric field 
is weaker. This is important for attempts to extract the 
dynamical critical exponents from experimental data^ 
This paper is organized as follows. In Sec.^we intro- 
duce our model and formulate an effective field theory. 
This theory is investigated in a Gaussian approximation 
in Sec. 11111 which yields a generalized diffusion equation 
for density fluctuations. In Sec. IIVI we derive a gener- 
alized nonlinear cr-model, which is investigated in a one- 
loop approximation. In Sec. we discuss our results 
and their connections to previous investigations. Some 
technical points are relegated to an appendix. 



II. MODEL, AND FIELD-THEORETIC 
FORMULATION 



We consider the Hamilton operator 



H — H^, 



V{x) 



Here 



A 

2m 



(2.1) 



(2.2) 



is the Hamilton operator for free particles in the presence 
of an electric field F, with A the Laplace operator and m 
the electron mass, and V{x) is a random potential. We 
use units such that h and the electron charge are equal 
to unity. V is characterized by a Gaussian distribution 
with zero mean and second moment 



{V{x)V{x'))i,. 



1 



6{x-x') 



(2.3) 



Here (. . .)dis denotes the disorder average, v is the den- 
sity of states per spin at the Fermi level, and r is the 
single-particle scattering time. Unless otherwise noted, 
we consider a two-dimensional system, d = 2. 

In our investigation we focus on configuration aver- 
ages of retarded and advanced Green functions and their 
products. These functions are solutions of the differential 
equations 



{±iuj + E-H) G"'^(a;, x'\E; oj) = S{x - x') 



(2.4) 



In equilibrium, F — 0, E = ^ is the chemical potential. 
For F ^ 0, the quantity E — F ■ x = fix the chemi- 
cal potential in a local equilibrium approximation. As a 
boundary condition, we require that the Green functions 
vanish at infinity. ^^'22, The solutions of Eq. (|2.4|l then are 
symmetric with respect to an interchange of x and x' , 



G'^-''^{x,x'\E;lu) = G^^'^{x',x\E;uj) 



(2.5a) 



This symmetry is a direct consequence of the time- 
reversal invariance, which is not broken by the electric 
fieldii 



The electric field does, however, break the translational 
invariance in real space. In the absence of the random 
potential, or after disorder averaging, a translation by a 
vector a leads to an energy change F • a. The config- 
uration averages G^'^ of the Green functions therefore 
satisfy the relationship 

G^^^{x + a,x' + a\E; uj) ^ G^'^(a;, x'\E -F ■a;uj) . 

(2.5b) 

The symmetry properties expressed by Eqs. (|2.5|l will be 
important later. 

We now consider a generating functional for our Green 
functions. Following Ref. we define 



where 



S"^"^ = j dx $(a;) {±iuj + E - H) <^{x) 



(2.6) 



(2.7) 



with $(a;) a real scalar field. The retarded and advanced 
Green functions are obtained according to the rule 



G^^^{x,x'\E;uj) 



(2.8) 

To calculate the configuration average we use the 
replica trick. To this end we consider 2n copies of the 
original generating functional (n for the retarded degrees 
of freedom, and n for the advanced ones), and take the 
limit n ^ at the end of the calculation. The replicated 
generating functional takes the form 



(z«z^)" = / n ^['^'"] 



Here a = — (n— 1) , . . . , n, A^, = —1 for a < and A^^ = 1 
for a > 0, uja = wAq, and ja — jA for a < and ja — jn 
for a > 0. Now we calculate the configuration average 
Z ~ ((Z^Z^)")dis, which generates the averaged Green 
functions G^''^. If one performs a Hubbard-Stratonovich 
transformation to decouple the resulting term quartic in 
<&, and then integrates out <I>, one finds 



Z[j = 0] = / D[Q] e 



(2.10) 



where 



^[Q] ^^Jdxtr {Qixjf - i tr InG^^ , (2.11) 



with 



{GQ^)aa'{x,x'\E;uj) = [{iuJa+ E - Hx)5aa' 



2t 



Qc.^.ix)]6{x-x') . (2.12) 



3 



It is easy to show that Gq , averaged with respect to the 
action A, and taken in the Hmit n — > 0, is diagonal in the 
repUca indices, and equal to {G^) for a > (a < 0). 
The Q-matrix fields in the Eqs. (|Tn|) and ^(TT^ are 
real-valued matrix fields which satisfy the relationship 



(2.13) 



The Eqs. (|2.11|) - H2.13|l define the effective action which 
forms the basis for our further investigations. 



III. TRANSPORT AND HEATING IN 
GAUSSIAN APPROXIMATION 

A. Saddle-point Green functions 

We now look for a saddle-point solution Q^^ of the 
effective action A[Q], defined by 



5A[Q] 



SQ 



(3.1) 



The saddle-point values of Q and Gq are diagonal in the 
replica indices, Q^^, = Saa' Qa, (Gqsp)qq' = ^gg' Gq, as 
are the exact expectation values. From Eqs. (|2.11l 
we find 



Qa{x) = — Ga{x,x\E;uj) 

TTV 

= — G^{0,0\E-F-x;io) , 

TTV 

= Qa^ip,) . (3.2) 

In the second equation we have used the generalized 
translational invariance property of the Green function. 



Eq. (|2.5b|l . and the notation in the third equation em- 
phasizes that the saddle-point field depends on x and 
the direction of F only via Hx = E — F ■ x. 

In order to obtain an equation for the Green function 
we use the generalized translation invariance property 
(|2.5bl) and the time reversal invariance. Due to these 
properties the Green function in the presence of the elec- 
tric field can be written in the formH 

G.{x,x'\E-u) = /(|pe--(— ') 

xg^{k\E-F ■{x + x')/2;uj) , (3.3) 

where ga{k\E]ijj) is a function that depends only 
quadratically on the field F . If we work to linear or- 
der in the field, the function ga{k\E\Lo) thus reduces to 
the equilibrium Green function. 



ga{k\E;u}) 



1 



E-ek + i^Q»{E) 



(3.4) 



where efc = fc^/ (2m). The only difference to the situation 
in the absence of the field is then that the quantity /i^:. 



which has replaced the chemical potential, goes to zero 
at the classical turning point, where F ■ x = E. We 
will nevertheless consider \ijLx\ large compared to 1/(2t) 
in calculating the integral in Eq. (|3.3|l . i.e., we consider 
only the region far from the classical turning point. In 
this approximation we obtain 



Qa{x) = sgnwg 



(3.5) 



In writing down Eq. I|3.5|l we have ignored the imagi- 
nary part of Qa, which only leads to a weakly field- 
dependent renormalization of the chemical potential. 
The saddle-point Green function obtained in this way 
agrees with that derived in Ref. 23 in a self-consistent 
Born-approximation. 

We note that the simple structure of the above saddle- 
point solution is a result of the constant density of states 
in our model. If the density of states were energy depen- 
dent, the saddle-point field would pick up an energy and 
field dependent contribution from the density of states, 
which would modify Eq. (|3.5|) . 



B. Gaussian fluctuations 

We now consider the Gaussian fiuctuations about our 
saddle-point solution. Substituting Q = Q^^ + 6Q into 
Eq. H2.11|l . we find for the Gaussian part of the action 



1 



87 



J dxdy'Y^ 6Qaa'{x) 



Six - y) 



2'KVT 



Ga'{y, x\E; uj) Ga{x, y\E; uj) 



<SQa'a(y) . (3.6) 



The properties of the Gaussian propagators depend on 
whether both Green functions in Eq. (|3.6|) are retarded 
or advanced, or whether one is retarded and the other 
advanced. In the former case, the propagator is massive, 
in the latter, soft. Setting ui — fl/2, the soft propagator 
P satisfies the equation 

J dy Fix, y\E; n) F(y, x'\E; n) = S{x - x') , (3.7a) 



where 



Tr{x,y\E;n) = 5{x-y)- 



1 



2'KVT 

xG^,<o{y,x\E;n/2) , 



G„>o{x,y\E;n/2) 
(3.7b) 



is the corresponding soft vertex. P describes the relax- 
ation of the particle number density.^ 

From Eqs. (|2.5() we obtain the following properties of 
the vertex function F, 



r{x,x'\E;n) = r{x',x\E;n) 



(3.8a) 



and 



F(a; + a,x' + a\E; ft) = F(a;, x'\E - F ■ a^fl) . (3.8b) 



4 



These relations are important for a gradient expansion of 
the vertex function. With their help, Eqs. H3.7|l can be 
written, to second order in the gradient operator and in 
the classically accessible region, 

{fl-V ■ D{fi^)V)P{x,x'\E;n) = S{x- x') . (3.9) 

Here 

D{nx) = IJ-xT/m (3.10) 

is the electron diffusion coefficient in a local equilibrium 
approximation. The details of the derivation are given 
in Appendix IXI We see that, inside the classically acces- 
sible region, the soft propagator is governed by a gener- 
alized diffusion equation, as one would expect from the 
fact that physically, P describes the relaxation of den- 
sity fluctuations. Outside the classically accessible re- 
gion, V{x,x'\E\Q) Ri 6{x — x'). Accordingly, only the 
modes inside the classically accessible region are general- 
ized diffusion modes, again in agreement with what one 
would expect. Equation l|3.9|) is very similar to the equa- 
tion derived in R.ef. ITtI bv kinetic theory methods. 

For later reference we note that, in a schematic nota- 
tion, the Gaussian propagator P has the form 

P =l/{n- DV^ + {D/E)F -V) . (3.11) 

Accordingly, there is a frequency scale 

n* = DF^/E^ , (3.12) 

or a time scale, t* = 1/^2*, that separates diffusive be- 
havior at f2 3> ri* from field-dominated drift behavior at 

n < n*. 



In deriving the generalized diffusion equation we have 
taken into account only terms linear in F. An estimate 
of the terms quadratic in F shows that they are small 
compared to the ones kept provided that \F\£/iJ. ^ 1, 
with £ the mean-free path. 



C. Density relaxation 

The generalized diffusion equation, Eq. H3.9() . differs 
from the ordinary diffusion equation by the real-space 
dependence of D, which produces a term linear in the gra- 
dient. This introduces a new singularity into the differ- 
ential equation. As a result, the solution of this equation 
for long times differs strongly from the ordinary Gaus- 
sian one. Let the field point in x-direction, F — (F, 0), 
put X = (x, y) and x' = (x' , 0)r^ and perform a Fourier 
transform with respect to the direction transverse to the 
field, 

P{x,k;x'\E;n) = J dy e''"' P{x,x'\E;n) . (3.13) 
We then find that the solution of Eq. (|3.9|l takes the form 



P{x,k;x'\E;n)^r{tx,,\k\;f^x'\E;n) , (3.14a) 
where 



v{^l,p■,^x'\E■,n) 



D'F 



-p{p+li')/F) 



e{n~fi')M -{l + n/D'Fp),l,2p^l'/F] U -{l + n/D'Fp),l,2pfi/F 



+e(/i' -^)u(-il + n/D'Fp), 1, 2p^i'/F Af -(1 + n/D'Fp), 1, 2pfi/F 



T[-il + n/D'Fp) 



(3.14b) 



Her e D' = D{E)/E = r/m. The functions U and M in 
Eq. (|3.14b|) are the confluent hypergeometric functionsi^i 
In order to obtain this solution we have required that P 
vanishes at infinity in the classically accessible region and 
that the probability current vanishes at the turning point, 
so that the particles can not penetrate into the classically 
forbidden region. 

To understand the nature of this solution we recall 
that P describes the relaxation of the particle number 
density. Let us consider, at time i = 0, an ensemble of 
particles with energy E = /i, located at x = 0, and with 



a homogeneous density in the direction perpendicular to 
the electric field. Accordingly, the initial number density 
is given by 

no{x,E)^'^5{E- ^i)5{x) , (3.15) 

where L is the linear dimension of the system in y- 
direction. The evolution of this initial density is governed 
by the propagator P for fc = 0. In this limit V takes the 
form 
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v{^JL,Q■,^Ji'\E■^) 



D'F 



1 4r2/i 

D'F^ 



1 4r2/i' 



(3.16) 



Here Iq and i^To are the modified Bessel functions. Tlie 
inverse Laplace transformation yields 

(3.17) 

Equation H3.17|l explicitly displays the characteristic time 

t* = E/D'F^ , (3.18) 

which was apparent already in Eq. H3.9I) (see Eqs. 1)3.111 
13.121) ). and which serves as the boundary between the 
short-time and the long-time behavior. In order to deter- 
mine the structure of the generalized diffusion propagator 
for t <t: t*, we use Eq. H3.14a|l . For t -^t* the spread of 
the initial (5-package in a;-direction is small, and we can 
use Fx/E^Fx'/E <C 1, in addition to the asymptotic 
expansion of /q for large arguments. If we expand the 
exponent with respect to a; — a;', we find 

P{x,Q]x'\E;t) ^ ^ 



^-{x-x'Y /iD{E) t 



y^4:7rD{E)t 

. ^^-^^^ 

As expected, the relaxation of the initial density per- 
turbation is diffusive in this time regime. Note that in 
writing Eq. H3.19|l we have ignored the first moment of 
the generalized diffusion propagator, which is nonzero. 
While the first moment is crucial for calculating currents, 
it is irrelevant for our current discussion. 

In the opposite limit, t ^ t*, the behavior is very 
different. In this case the width of the particle packet 
becomes very large and the asymmetry of the particle 
packet, which is small initially, is getting considerable. 
The Bessel function in Eq. (|3.17|) approaches unity for 
asymptotically long times, so that 

^(--0;x'|i?,t) = -^ exp(-^^^) , (3.20) 

for ^HxfJ^x' /E>'F^t <C 1. For very large \x\, which satisfy 
the requirement ^/JT^JT^/D'F'^t 1, we obtain 



Pix,0;x'\E, t) = 



1 



exp(- 



/fJ.x 



-) 



D'Ft D'FH 

(3.21) 

These results show that in this regime the dynamics is 
no longer diffusive. 

If, instead of the distribution function H3.15|l . we con- 
sider a delta pulse in both the longitudinal and trans- 
verse directions, we have to investigate the function P 
for nonzero values of k. We have not been able to do so 
exactly, and have resorted to a WKB approximation in- 
stead. The results obtained in this way are in qualitative 
agreement with the case analyzed above. 



D. Transport and heating 

So far we have considered the relaxation of density per- 
turbations. We now turn to the questions of transport, 
and Joule heating. To this end we consider, instead of 
Eq. H3.15|l . an initial density distribution 



nQ{x,E) 



2'kN 
V 



5(E- ii-F -x) 



(3.22) 



Here TV is the total particle number, and V is the sys- 
tem volume. Notice that Eq. 1)3.22(1 describes a uniform 
number density in the classically accessible region, since 

/dE 
— no{x,E)^e{fi,)N/V . (3.23) 

To see how such a density distribution evolves, we re- 
alize that V, Eq. H3.17|l . gives the probability to find a 
particle at energy fi' at time t, if it had the energy fi at 
time t = 0. In the limit of long times, t ^ t*, the mean 
particle energy therefore increases linearly with time, 

= fi + D'F^t . (3.24) 

Notice that P and V were normalized with respect to 
an integration over x. Changing the integration variable 
to fi results in the additional factor 1/_F in Eq. H3.24|l . 
Equation (|3.24|) shows that the energy fed into the system 
increases the kinetic energy of the charge carriers. 

According to the generalized diffusion equation, the 
current density distribution is given by 



j{x,E\n) = -D{fi,)Wn{x,E\n) 



(3.25) 



where n(x, E\Q) is to be calculated from no by means of 
the propagator P. The volume averaged current density 
takes the form 

i(f^) = l J ^dxj{x,E\n) . (3.26) 

For our spatially uniform charge density distribution, 
n{x,E\Q) = n(fix\^), and we obtain 

d/i dD(fi) 



(3.27) 



2tt dfi 

In our Gaussian approximation, the derivative of the dif- 
fusion constant is independent of /it, and thus can be 
taken out of the integral. We finally obtain 

. . rN F , , 

The Gaussian theory thus yields an Ohmic current that 
leads to Joule heating. 
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IV. NONLINEAR a-MODELL 
A. Effective action 

The derivation of a matrix field theory in Sec. Illll has 
proceeded in analogy to the case without an electric field, 
and the result was structurally very similar to the latter. 
In particular, matrix elements that correspond to prod- 
ucts of retarded and advance degrees of freedom are soft, 
while those corresponding to products of two advanced or 
two retarded degrees of freedom are massive. The chief 
difference is that, in the presence of an electric field, the 
soft modes in Gaussian approximation are not diffusive, 
but rather obey the more complicated differential equa- 
tion 

It is obvious from these observations that one can de- 
rive an effective field theory for the soft modes by repeat- 
ing the procedure that leads to a nonlinear cj-model in 
the zero-field casei^^i^. As expected, the only change is 
that the Laplace operator in the cr-model is replaced by 
the differential operator from Eq. (|3.9|l . Wc thus find 



n dxtr (AQ{x) 



Here A is a diagonal matrix whose elements have been 
given after Eq. H2.9|) . The matrix elements of Q are el- 
ements of the homogeneous space 0{n,n)/0{n) x 0{n), 
and Q is subject to the constraints 



trQix) ^ , 
Q^{x) = 1 . 
Q^{x) = Ag(a;)A 



(4.2a) 
(4.2b) 
(4.2c) 



These can be incorporated in a parameterization in terms 
oi n X n matrices q. 



(4.3) 



Roughly speaking, the field Q contains the soft parts of 
the field Q of the previous section. 

An expansion of the action, Eq. H4.1|l , to quartic order 
in q reads 



(4.4a) 



(4.1) with 



^(2) 
^(4) 



TTI/ 

32 



^ J dx^qaa'ix) (n - V ■ D{fix)'V^ qaa'{x) , 

Oca' 

J dx ^qaia2{x) qaaaiix) (j^ - V ■ D {fl^)'V^ qasaix) qa^aA^) 



Qa4a3 {x) qaicti {x) 

r 



(4.4b) 



(4.4c) 



B. One- loop theory for the diffusivity 

The expansion of the effective action in powers of q 
allows for a systematic loop expansion. To one-loop order 
we find the following result for the diffusivity D, 



1 - — P{x,x\E;n) 



, (4.5) 



which has the same structure as in the absence of the 
electric field. The only difference is that P{x, x\E, il) is 
calculated from Eq. (j3.14b|) , according to the relationship 



p{x,x\E,n)^ [ ^v{^i,,p;^i.\E,n) , (4.6) 

Jo 

where A is an ultraviolet cutoff. In the absence of the 
electric field this integral is infrared divergent if = 0. 



In the presence of the electric field we obtain, for = 0, 
and in the limit Xfj,x/F ^ 1, 



(4.7) 

The static one-loop diffusivity is thus finite in two- 
dimensions, indicating that the electric field destroys the 
mechanism that gives rise to localization in d = 2. This 
was already obvious from Eq. I|3.11|l . which is less infrared 
divergent than a diffusion propagator. 

For n > 1/r, with t* from Eq. (jHUHJ, the function 
P reduces to the conventional diffusion propagator, as 
discussed in Sec. IIIII Consequently, in this regime the 
corrections to the bare diffusion coefficient take the same 
form as in the conventional weak-localization theory in 
the absence of the field. 
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C. Scaling analysis 

We now perform a scaling analysis of the generalized 
nonlinear cr-model. Our procedure is analogous to the 
one in Ref. IsS To this end, we deviate from our restric- 
tion to d = 2 and consider the model in c? > 2, ignoring 
the complications that arise from the chemical potential 
dependence of the density of states in d > 2. We first 
write the effective action in a schematic form that leaves 
out everything that is not necessary for power-counting 
purposes, 

A^a = J dx V^q^ + H J dx n 

+ 1- f dx xV^q^ + ^ [ dxV^q'' + ... (4.8) 

Here iJ cx j/, I/G2 oc i/r, etc. We now assign a scale 
dimension [L] ~ —1 to lengths L. For F = 0, this action 
contains two fixed points, namely, a stable Gaussian one 
that describes the diffusive phase, and a critical one that 
describes the Anderson transition. 



1. Gaussian fixed point 

The Gaussian fixed point that describes diffusion in 
the absence of an electric field, F = 0, is characterized 
by a scale dimension of the field q equal to 

Mdiff = (d - 2)/2 . (4.9a) 

Frequencies must scale like wavenumbers squared in a 
diffusive phase, so we also require 

[J^]diff = 2 . (4.9b) 

G2 and H are then dimensionless, 

[G2]diff = [-ff]diff = , (4.9c) 

and all higher order terms are irrelevant, with the least 
irrelevant couplings having a scale dimension 

Mdiff = -(d-2) , (4.9d) 

where u denotes a generic irrelevant operator. I/G4 is 
an example of such a least irrelevant operator. Adding 
the electric field, we see that F is relevant with respect 
to the diffusive fixed point, 

[F]diff=l . (4.10) 

The frequency dependent diffusivity, whose bare value is 
D = I/G2H, therefore obeys a scaling law 

D{n,F,u) ^ D{nP,Fb,ub-^'^-^^) . (4.11) 

At = 0, and for small F, we conclude that the diffu- 
sivity has the structure 

D{F) cx const, -t- F'^~^ , (4.12) 



and in d = 2 one expects a logarithmic dependence of D 
on F. This is in agreement with the explicit perturbative 
result in the previous subsection. 

For 7^ 0, Z? is a function of F'^ /ft, in agreement with 
the explicit Gaussian theory in Sec. IIIII There are two 
scaling regimes. For small il, the scaling of D is gov- 
erned by the electric field, and for large fl, the scaling is 
governed by the frequency. The crossover between these 
two scaling regimes is at a frequency J7 = fi*, Eq. H3.12|l . 
In an experiment, the frequency O is effectively replaced 
by 1/t0, where is the phase relaxation time*^ The 
weak-localization physics can therefore only be observed 
if the temperature is large compared to a crossover tem- 
perature T*, which is given by 

1/t4,(T*) = DF^/E^ . (4.13) 

We will further discuss this condition in Sec. IVl below. 



2. Critical fixed point 

We now turn to the critical fixed point that describes 
an Anderson transition at F = 0. Here we choose the 
field q to be dimensionless, and the scale dimension of 
the frequency to be d, 

[q],^0 , [n], = d . (4.14) 

The bare scale dimension of G2 is then [G2]c = [GaIc — 
. . . = [G]c = 2 — d = — e. An explicit renormalization- 
group calculation shows that H is not renormalized, while 
the renormalized counterpart of G, g, has a fixed point 
value g* = 0{e). The deviations of g from g* constitute 
the relevant operator at the critical fixed point, whose 
scale dimension determines the correlation lenght expo- 
nent v. To one-loop order— 

i^ = l/e + 0(l) . (4.15) 

The scale dimension of F is also given by a loop expan- 
sion, but the leading term can again be determined just 
by power counting, 

[F], = l + 0(e) . (4.16) 

The electric field is thus a relevant operator with respect 
to the critical Anderson fixed point. In particular, an 
arbitrarily weak field F destroys the usual localization in 
d = 2. 

The above discussion shows only that the F — fixed 
point is unstable, and does not tell what happens instead. 
The perturbative result, Eq. 14.7|l . suggests that there is 
a metallic phase in d = 2 for F 7^ 0. We note that the 
result, Eq. (|4.16|) . is different from the popular scaling 
argument which assumes that F^, with ^ the correlation 
length, represents the critical energy or frequency scale, 
which yields [F] = d -I- lii^ We will come back to this 
discrepancy in Sec. IVl below. 
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V. DISCUSSION 

In summary, we have used field-theoretic methods to 
investigate the impact of an electric field on the lo- 
calization of noninteracting electrons, mostly in two- 
dimensions. We have found that there is a character- 
istic temperature T* that separates a regime where the 
physics is dominated by the field from one where it is not. 
For T < T* the physics is dominated by the electric field, 
which directly affects the structure of the localization cor- 
rections. In this regime the density relaxation is strongly 
non-diffusive, and the usual weak-localization corrections 
to observables are replaced by logarithmic dependences 
on the field. The latter have the same structure as those 
derived in Refs.l20landl23. Our treatment shows that the 
applicability of these results is restricted to T < T* . For 
T > T* the electric field does not significantly affect the 
diffusion of a particle packet. Consequently, the weak- 
localization corrections to the diffusion coefficient in this 
regime are the same as in equilibrium, and independent 
of the electric field. In this regime the approaches of Refs. 
El and m are well founded. 

Let us estimate the value of T* for parameter val- 
ues that are representative of a typical weak-localization 
experiment,'^ With £> « 14 cm^/s, F « 1.6 x 
lO^^eV/cm, and E « 0.6 meV, one has 1/t^{T*) « 
lO^Hz. At low temperatures, is dominated by the 
electron-electron interaction and inversely proportional 
to the temperature,-- = c/T. For the data of Ref. fssi 
we obtain c « lO^K-^s^^ This yields T* « 10"^ K. We 
conclude that the crossover between the field-dominated 
regime and the usual weak-localization regime occurs at 
unobservably low temperatures. This explains why no 
field-dominated scaling corrections were observed in the 
experiments of Refs. andll^i and it shows that these 
observations are not at odds with the notion that an ar- 
bitrarily weak electric field does indeed destroy the weak 
localization. 

We now come back to the scale dimension of the field 
F with respect to the critical fixed point that describes 
the metal-insulator transition. As mentioned in Sec. IIVI 
a scaling argument given by Sondhi et alii^ assumes that 
F£^ scales like the critical energy scale. For the scale di- 
mension of F this yields [F]c = z + 1, with z the dynam- 
ical critical exponent. Schematically, this corresponds to 
a critical propagator (in the case of an Anderson transi- 
tion) 



n + D{v)v^ + F/v ■ ^ ■ ^ 

This is not consistent with perturbation theory, however, 
at least not for the model studied here. From Eq. (|3.11(l 
we see that the actual critical propagator has the struc- 
ture 



n D(V)V2 -f Fi:'(v)v ■ ' ■ ' 

The difference between these two expressions accounts for 



the difference between the present result [F]c = 1, and 
that of Ref. IS. The reason for the additional gradient 
squared in Eq. 1)5.2(1 compared to Eq. ((5.1|l is the Ward 
identity"* that is closely related to particle number con- 
servation. Equation 1(5. 1() implies that the electric field 
breaks particle number conservation, which it does not. 
The same point can be made at the level of a fermionic 
action,^ which can be seen as follows. The F ■ x term in 
the Hamiltonian, Eq. (|2.2|l . corresponds to a term 

Sf = j '^^ dr {F ■ x) n{x.T) (5.3) 

in the action, with n the electron number density field. 
The latter corresponds to n = tr Q in the Q-TnaXriyi for- 
mulation of the field theory, and in the nonlinear a-model 
this corresponds to tr Q = 0. This shows that the leading 
coupling of F to the electronic soft modes vanishes for 
symmetry reasons. The leading nonvanishing term car- 
ries an additional gradient squared. This is the reason 
why one cannot obtain the nonlinear cr-model for elec- 
trons in an electric field by simply replacing the coupling 
term in the fermionic action by its Q-field counterpart. 
In contrast, an external magnetic field does break a sym- 
metry (viz., spin rotational invariance) and gives some 
soft modes (viz., the spin-diffusons) a mass, and replac- 
ing the spin density in the Zeeman term by its corre- 
sponding trace over a Q-field gives the correct answer. 
We conclude that, in general, one cannot extract values 
for the dynamical exponent z from experimental data by 
assuming [F] — z as was done in Ref. |28[ 

At this point we would like to note that our model 
does not take into account the electron-electron interac- 
tion and the resulting inelastic collisions. Therefore, our 
results are only valid for sample sizes smaller than the en- 
ergy relaxation length. If inelastic collisions were taken 
into account, the electron system should be described 
by a distribution function with an effective temperature. 
If the energy-transfer rate within the electronic subsys- 
tem is larger than the one between the electron and the 
phonon system, the effective temperature Toff would de- 
pend also on the strength of the applied electric field. 
Since in this case the temperature T in the in the phase 
relaxation time t^(T) would be replaced by Tes, the field 
dependence of the effective temperature would also be 
reflected in the logarithmic corrections to the conduc- 
tivity, as pointed out in Ref. ISTl This can be mistaken 
for a direct impact of the electric field on the localiza- 
tion corrections in a regime where there is actually none. 
However, these effects depend on the ratio between the 
energy-transfer rate in the electronic subsystem and the 
cooling rate. They are therefore not universal, but rather 
dependend on the experimental conditions. 

We finally discuss the relation of our theory to some 
previous work in the literature. We have already noted 
that our equation for the crucial Gaussian propagator, 
Eq. H3.9|) . is very similar to the one obtained by Kirk- 
patrick with different methods^ Indeed, our theory is 
in many respects a field-theoretic version of his kinetic 
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theory. A different differential equation for the propaga- 
tor was used in Refs. i2L22.23i. which found that a critical 
field strength Fc is necessary to delocalize the states at 
the Fermi energy. Instead of Eq. (|3.9|) , these authors used 



{n ~ D{E)V^ + D'{E) F ■ V) P{x, x'\E; Q) = S{x-x'), 

(5.4) 

with D'{E) = dD{E)/dE. The use of this equation was 
based on the notion that there is a rapid mechanism for 
cooling, so that heating processes can be ignored. The 
resulting equation violates the properties expressed by 
Eqs. H3.8|l . which are based on time reversal symmetry 
and generalized translational invariance. Indeed, if there 
is a rapid mechanism for cooling that validates Eq. H5.4|l . 
then the electrons must obviously have already experi- 
enced inelastic collisions before the time at which Eq. 
H5.4|l becomes valid. They are thus already in the diffu- 
sive regime, t > t* , where the electric field no longer pro- 
vides the leading effect. In the opposite regime, t < t*, 
Eq. H5.4|) is not valid, and Eq. H3.9|) must be used instead. 
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APPENDIX A: DERIVATION OF THE 
GENERALIZED DIFFUSION EQUATION 

In order to obtain Eq. (|3.9|) we have to expand F about 
the (5-function. To this end we need to calculate the mo- 
ments of F. We first consider the zeroth moment, which 
is given in terms of 

Mo = j dvGo,{x,y\E-oj)G-^{y.x\E;oj) , (Al) 

with Lo — and a > 0. In terms of the function g^n 
Eq. H3.4|l . we have 

Mo = I dy I j g^{k\E~F ■{x-y/2);uj) 
x.g_„(p|i?-F.(x-y/2);c.)e-(''-P>'' , (A2) 



where = /dfc/(27r)2. The terms F ■ y/2 in Eq. (|X2|l 

lead to corrections of 0{F^). If we omit them we obtain 

Jk 

de 1 



+ jpsgn (w) iuj-ek+ jT^gn (-w) 



(A3) 

In order to evaluate this integral we use the approxima- 
tion already discussed in Sec. IIIII li fi^ > we replace 
fix by cxD, and if fix < we replace fix hy —oo. In the 
hydrodynamic limit we then obtain 



A/o = 2TTVTQ{fj.x) (1 - ^t) 



The first moment Mi is defined as 



(A4) 



Ml = J dy {x- y) G^{x, y\E- lo) G-c.{y, x\E; uj) 



In the same approximation as above we find 



(A5) 



Ml 



F d 



4 dE 

F-^M2 . 

dE 



y{k + K\fix; uj )g-.a{k\fix;uj) 



(A6) 



Here M2 is the second moment. It is defined as 

M2 = \j dy {y - xf Go.{x,y\E;u)G-^{y,x\E-uj). 

(A7) 

Here we ignore the fact that anisotropic terms can arise 
due to the electric field, since any such terms are of 
0{F^). Again using the same approximations as above, 
we obtain in the hydrodynamic limit 



Mo = 



2m 



1 



de , , S — = 2TTVTD{fix)T9{fix). 

(A8) 

Collecting our results, we now have the following expres- 
sion for F in the hydrodynamic regime, and for x in the 
classically accessible region, 

T{x, y) ^5{x-y)[n-V ■ D{fix)V] . (A9) 

From this, Eq. (|3.9|) follows immediately. 
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